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We use QCD sum rules to study the possible existence of a D S D* + D^D molecule with the 
quantum number J p = 1 + . We consider the contributions of condensates up to dimension eight 
and work at leading order in a s . We obtain uid b d* = (3.96±0.10) GeV around 100 MeV above the 
mass of the meson X(3872). The proposed state is a natural generalized state to the strangeness 
sector of the X(3872), which was also found to be consistent with a multiquark state from a previous 
QCD sum rule analysis. 
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The recent appearance of new resonances observed by BaBar, BELLE, CLEO and FOCUS collabora- 
O i' tions, and partly confirmed by the Fermilab experiments CDF and DO, sheds new light on the spectroscopy 

of charmonium states. Among these new mesons, some have their masses very close to the meson-meson 
threshold like the X(3872) [l| and the Z + (4430) Q. Of special importance is the appearance of the 
Z + (4430), which decays into ip' and n + and, therefore, can not be described as ordinary cc meson. 
Its nature is completely open, but an intriguing possibility is the interpretation as tetraquark state or 
molecular state 0, 0, 0- If these mesons are really molecular states, then many other molecules should 
exist. A systematic study of these molecular states and their experimental observation would confirm its 
structure and provide a new testing ground for QCD within multiquark configuration. In this context 
a natural extension would be to probe the strangeness sector. In particular, in analogy with the meson 
A(3872), a D S D* molecule with J p = 1 + would decay into J/ipK* — > J/ipKn and, therefore, could be 
easily reconstructed. Here we use the QCD sum rules [1, 0, H| to predict the mass of such a molecular 
jr;* ■ state. 

This calculation is of particular importance for new upcoming experiments which can investigate with 
much higher precision the charmonium energy regime, like the PANDA experiment at the antiproton- 
proton facility at FAIR, or a possible Super-B factory experiment. Especially PANDA can do a careful 
, scan of the various thresholds being present, in addition to precisely going through the exact form of the 

5_i ■ resonance curve. 

03 ' In a previous calculation, the QCDSR approach was used to study the A (3872) considered as a diquark- 

antidiquark state 0, and the Z + (4430) meson, considered as a D*D\ molecular state @. In both cases 
a very good agreement with the experimental mass was obtained. 

Considering a D*D S molecule with J p = 1 + , a possible current describing such state is given by: 

i 

3n= —^[{Sal5Ca){cblf,d b ) - (s a "f^C a ){c b ~f 5 d b )} , (1) 

where a and b are color indices. We have considered the anti-symmetrical state D*~ Df — D~ D s * + to keep 
a closer relation with the A(3872) meson. Considering the A(3872) as a D*D molecule, the combination 
D*°£) — D°D*° has J PC = 1 ++ as the A(3872) meson. Of course the symmetrical combination: 
D*~ Df + D~ D s * + would provide exactly the same mass, within our sum rule approach. 
The sum rule is constructed from the two-point correlation function: 

IW?) = ijd^x e^(0\T[j,(x)ji(0)}\0) = -Mq 2 )(g» v - ^) + n (g 2 )^. (2) 
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Since the axial vector current is not conserved, the two functions, III and IIo, appearing in Eq. @ are 
independent. They have respectively the quantum numbers of the spin 1 and mesons. Therefore, we 
choose to work with the Lorentz structure g^, since it gets contributions only from the 1 + state. 

On the OPE side, we work at leading order in a s in the operators and consider the contributions from 
condensates up to dimension eight. The correlation function, III, in the OPE side can be written as a 
dispersion relation: 



n? PE tf) = / ds 

J 4m? 



oo JDPE 



4m? s - 1 



(3) 



where p OPE {s) is given by the imaginary part of the correlation function: irp OPE (s) = Impi (s)]. We 
get: 

p OPE {s) = p per \s)+p m '(s)+p {Sg \s)+p m ^ q \s)+p {G2 ^s)+p mix (s)+p^ , 



with 
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where the integration limits are given by a r< 



(1 - Vl - 4m2/s)/2, a max = (1 + y/l - 4m2/ s )/2, 



/?min = ar7i^/(sa — m 2 ), and we have used (qga.Gq) = m^qq), (sGa.GS) = mffi. 

It is very interesting to notice that the current in Eq. |T]) has a OPE behavior very similar to the 
scalar-diquark axial-antidiquark current used for the A (3872) meson in ref. 0]. 

In the phenomenological side, we write a dispersion relation to the correlation function in Eq. (U): 



phen 



ds 



(G) 



where p phen (s) is the spectral density and the dots represent subtraction terms. The spectral density is 
described, as usual, as a single sharp pole representing the lowest resonance plus a smooth continuum 
representing higher mass states: 



P phen {s) = \ 2 6(s-m 2 DsD ,) 
where A gives the coupling of the current to the meson D S D* 

(0\j\D s D*) = A. 



(7) 



(8) 



For simplicity, it is assumed that the continuum contribution to the spectral density, p cont (s) in Eq. (O, 
vanishes bellow a certain continuum threshold s . Above this threshold, it is assumed to be given by the 
result obtained with the OPE. Therefore, one uses the ansatz [lOl 



n OPE 



(s)9(s - s ) , 



(9) 



After making a Borel transform to both sides of the sum rule, and transferring the continuum contribu- 
tion to the OPE side, the sum rules for the pseudoscalar meson Z + , up to dimension-eight condensates, 
can be written as: 



AV 



ds e' s l M p OPE {s) 



(10) 



To extract the mass rriD s D* we take the derivative of Eq. (JTUJ) with respect to 1/M 2 , and divide the 
result by Eq. (fTOjl . 
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FIG. 1: The OPE convergence in the region 1.9 < M 2 < 3.0 GeV 2 for ^ = 4.5 GeV. We start with the relative 
perturbative contribution (long-dashed line), and each other line represents the relative contribution after adding 
of one extra condensate in the expansion: +m s (dashed line), + (qq) (dotted line), + m a {qq) (dot-dashed line), 
+ (g 2 G 2 ) (line with circles), + rnx,{qq) (line with squares), + {qq} 2 (line with diamonds), + m s mo(gg) (line with 
triangles up), + m s (qq) 2 (line with triangles down), + m^qq) 2 (solid line). 

The values used for the quark masses and condensates are d, [ll|: m c {m c ) = (1.23 ± 0.05) GeV, 



(qq) = -(0.23 ± 0.03) 3 GeV , 8 = 0.8(qq), 



0.8 GeV , (g 2 G 2 ) = 0.88 GeV 4 
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We evaluate the sum rules in the Borel range 1.9 < M 2 < 3.5 GeV 2 , and in the sq range 4.4 < y/J^ < 4.6 
GeV. 

From Fig. Q] we see that we obtain a reasonable OPE convergence for M 2 > 1.9 GeV 2 . Therefore, 
we fix the lower value of M 2 in the sum rule window as M 2 nin = 1.9 GeV 2 . We notice that the OPE 
convergence in this case is similar to the OPE convergence for the X(3872) meson Q, and not so good 
as the OPE convergence for the Z + (4430) meson Q. 

0.9 

0.8 
0.7 

g 

E 0.6 

I 0.5 
o 
o 

» 0.4 
o 

Q_ 

0.3 
0.2 
0.1 

1.8 2 2.2 2.4 2.6 2.8 3 3.2 

M 2 (GeV 2 ) 

FIG. 2: The dashed line shows the relative pole contribution (the pole contribution divided by the total, pole 
plus continuum, contribution) and the solid line shows the relative continuum contribution for s /sq = 4.5 GeV. 

To get an upper limit constraint for M 2 we impose that the QCD continuum contribution should 
be smaller than the pole contribution. The comparison between pole and continuum contributions for 
y/so = 4.5 GeV is shown in Fig. O From this figure we see that the pole contribution is bigger than the 
continuum for M 2 < 3.0 GeV 2 . The maximum value of M 2 for which this constraint is satisfied depends 
on the value of Sq. The same analysis for the other values of the continuum threshold gives M 2 < 2.8 
GeV 2 for ^/s^ — 4.4 GeV and M 2 < 3.2 GeV 2 for ^/s^ = 4.6 GeV. In our numerical analysis, we shall 
then consider the range of M 2 values from 1.9 GeV until the one allowed by the pole dominance criteria 
given above. 

In Fig. [31 we show the D S D* meson mass, for different values of y/so, in the relevant sum rule window, 
with the upper validity limits indicated. From this figure we see that the results are very stable as a 
function of M 2 . 

Using the Borel window, for each value of sq, to evaluate the mass of the D S D* meson and then varying 
the value of the continuum threshold in the range 4.4 < ^/so < 4.6 GeV, we get 

m D ,D' = (3.97 ±0.08) GeV, (11) 

around 100 MeV bigger than the mass of the V(3872) meson. The Borel curve for the mass is quite 
stable and has a minimum within the relevant Borel window. Such stable Borel curve strongly suggests 
that there is indeed a very well defined ground state. 

To check the dependence of our results with the value of the charm quark mass, we fix y^so = 4.5 GeV 
and vary the charm quark mass in the range m c = (1.23 ± 0.05) GeV. Using 1.9 < M 2 < 3.0 GeV 2 we 
get: rriD s D» = (3.96 ± 0.10) GeV, in agreement with the result in Eq. (Til"j) . Therefore, we conclude that 
the most important sources of uncertainty in our calculation are the values of the continuum thrcshod 
and the charm quark mass. 

In conclucion, we have presented a QCDSR analysis of the two-point function for a possible D S D*+D*D 
molecular state with J p = 1 + . This state would decay into J/ipK* — > J/^Kir and, therefore, could be 
easily reconstructed. Our finding strongly suggests the possibility of the existence of such molecular 
resonance, whose structure is similar to the X(3872), with a mass about 100 MeV above that of X(3872). 

The mass of our proposed state is also close to the newly discovered X(3940) [13] and Z(3930) (l3j . 
While our proposed state has stangeness, the two states do not. Moreover, while V(3940) was found in 
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FIG. 3: The D S D* meson mass as a function of the sum rule parameter (M 2 ) for y/so = 4.4 GeV (dashed line), 
y/Io = 4.5 GeV (solid line) and y/Io = 4.6 GeV (dot-dashed line). The crosses indicate the upper limit in the 
Borel region allowed by the dominance of the QCD pole contribution. 

the decay of D*D and is a candidate for r\ c (3 1 Sq), .£(3930) was found in DD and is consistent with x' C 2- 
Therefore, a comparative study with these particles can give us clues on the structures of A (3872) and 
our proposed state. 

As a final remark, one can notice that if we take m s = 0, one obtains the sum sule for the A(3872) 
meson, considered as a molecular state D°D*° — D*°D a . The results are very similar to the ones obtained 
in ref. @, but more stable as a function of the Borel mass. Using y/so — (4.4 ± 0.1) GeV one obtains 
n%x = (3.88 ± 0.06) GeV, in an excelent agreement with the experimental value. 
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